Typical entanglement in multi-qubit systems 
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Quantum entanglement and its paradoxical properties hold the key to an information process- 
ing revolution. Much attention has focused recently on the challenging problem of characterizing 
entanglement. Entanglement for a two qubit system is reasonably well understood, however, the na- 
ture and properties of multiple qubit systems are largely unexplored. Motivated by the importance 
of such systems in quantum computing, we show that typical pure states of A'^ qubits are highly 
entangled but have decreasing amounts of pairwise entanglement (measured using the Wootter's 
concurrence formula) as A'^ increases. Above six qubits very few states have any pairwise entangle- 
ment, and generally, for a typical pure state of A'' qubits there is a sharp cut-off where its subsystems 
of size m become PPT (positive partial transpose i.e., separable or only bound entangled) around 
A'' > 2m + 3, based on numerical analysis up to A' = 13. 

PACS numbers: 03.67.-a, 03.65.Ud 

Quantum entanglement is a key prediction of quantum 
mechanics and is generally thought to be one of the cru- 
cial resources required in quantum information process- 
ing. Known quantum information applications include 



vice versa permuting A and B). For a mixed state p of 
two qubits, the concurrence C — r^/^ is given by 



C = max(Ai — A2 — A3 — A4, 0), 



(1) 



quantum computation ||^, 
PI , quantum cryptography 




quantum communication 
and quantum tele- 



portation |9|, |T^. The properties of entangled states 
potentially used in these applications are still poorly un- 
derstood. General two qubit entangled states have been 
well characterized, and a number of analytical measures 
of entanglement are known However, for N qubit 

systems (with iV > 2), few such measures can be cal- 
culated even for pure states. Despite these difficulties, 
arrays of qubits have been the focus of recent attention 
p3| , p^ , |l5|, |l^, e.g., Raussendorf and Briegel |17 pro- 
pose an array of highly entangled qubits for quantum 
computation. 

In this letter we consider the average, or typical proper- 
ties of N qubit pure states. A subsystem of an entangled 
pure state is in general a mixed state, thus we also con- 
sider the entanglement of mixed states derived from these 
pure states. Currently, a variety of measures are known 
for quantifying the degree of entanglement, including the 
entanglement of distillation the relative entropy of 
entanglement jl9j , the entanglement of formation |jri|, 
and the negativity [n2|. Using |0) and |1) for the two 
eigenstates of a spin^^ system or equivalent that encodes 
a qubit, and |00) as a shorthand for |0)a|0)s, a state of 
two qubits A and B, the four Bell states (|00) ± |ll))/\/2 
and (|01) ± |10))/\/2 have the maximum possible entan- 
glement. Most entanglement measures assign a value of 
1 to a Bell state, and to all separable states. 

For our purposes, it is convenient to use the tangle 
T (squared concurrence) as our entanglement measure 
when considering two or three qubits. The tangle is an 
entanglement monotone from which the entanglement of 
formation can be calculated 1 11 , ^ . For a pure state of 
two qubits, r = 4det|pyi|, where pA is the reduced density 
matrix obtained when qubit B has been traced over (or 



where the Ai are the square roots of the eigenvalues 
oi pp = p <7y ® Uy p*ay ® , AuA dcuotcs the 
complex conjugation of p in the computational basis 
{|00), |01), |10), 1 11)}. T he entanglement of formation 
Ef = ^(^ + — C*^), where h{x) is the binary en- 
tropy function, h(x) = — a;log2(a;) — (1 — a;) log2(l — xV 
A general pure state of two qubits can be written ||21| 



\^^) = e''^" cos 6*0 1 00) 

+ e*'^^ sin 6*0 cos 6*1101) 

+ e 

+ e 



'/^^ sin 6*0 sin6'i cos 6*21 10) 
^^"^ sin 00 sin6'i sin6l2|ll), 



(2) 



where ^re chosen uniformly according to 

(2^)-4d(sin0o)'d(sin0i)4d(sin02)'dxodxidx2dx3, (3) 

the Haar measure, with < < 27r and < 6*^ < 
7r/2. (We include an extra overall random phase, e*-*^", to 
maintain consistency with SU(7i).) 

Calculating t — ^A&\\pa\ using Eq. (||) and integrat- 
ing over the Haar measure, Eq. (^), gives (t) = 2/5. A 
randomly selected pure state of two qubits might thus be 
expected to have 0.4 tangle units of entanglement, and 
we have already noted that states exist with the max- 
imum (= 1) and minimum (= 0) amounts of entangle- 
ment. More informative is the distribution, the density 
of states (over the Haar measure) with a given value of 
the tangle. Calculated numerically by sampling 30 mil- 
lion random pure states, this distribution is shown in Fig. 
|l|. In tangle units, the distribution is broad, with many 
states having little or no entanglement compared to only 
a few with high or maximal entanglement. 

However, it must be emphasized that the shape of the 
distribution is heavily dependent on the choice of entan- 
glement measure. Compare in Fig. |l]the distributions of 



2 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

entanglement 

FIG. 1: Distribution of the tangle (*), concurrence (+) and 
entanglement of formation (■) for two qubit pure states sam- 
pled uniformly over the Haar measure. 



the concurrence C — r^/^ and the entanglement of forma- 
tion (first shown in |22]). The corresponding average val- 



ues (calculated numerically and in agreement with |22 ) 
are (C) = 0.5890 ± 0.0002 and {Ej) = 0.4808 ± 0.0002. 
Given the plethora of entanglement measures in current 
use, it is important to only compare like with like. Also, 
it is interesting to note that, though close, none of these 
average entanglement values are exactly equal to 1/2, 
which would be the nai've guess. 

For three qubits, it is possible to define the 3-tangle 
T3 for a pure state ||20| , giving a measure of the purely 
three-way entanglement in the system. A value for the 
tangle between each of the three possible pairs of qubits 
can also be calculated, t^B: ^ac and r^c? using Eq. (|^). 
These satisfy, 



TZ — TA- TAB — TAG 



(4) 



where ta = 4det|py!i| as before, except now both qubits B 
and C have been traced out to leave the partial density 
matrix pA- Equation (^) holds for any permutation of 
B and C. For the GHZ state (|000) -t- |111))/V2, which 
has the maximum possible 3-tangle, T3 = 1 and tab ~ 0, 
while for the W state (|001) + |010) + |100))/v%, T3 = 
and TAB = 4/9 for each pair, the maximum possible 
amount of pairwise tangle in a three qubit state |l5| . 

As for two qubits, the average values for ra and tab can 
be calculated. After analytically evaluating somewhat 
lengthy integrals, they are found to be (ra) = 1/3 and 
(tab) — 1/6- Since T3 and tab satisfy Eq. (^), one can 
consider the quantity (ra -I- tab + tag + tbc) =5/6, the 
average total entanglement (in tangle units) in a random 
pure state of three qubits. This gives an average entan- 
glement of 5/18 ~ 0.27 per qubit for three qubits, com- 
pared with 1/5 per qubit for two qubits. A GHZ state has 
a tangle of 1/3 per qubit and a W state 4/9. O'Connor 
and Wootters ||l^, |l^ considered rings and chains of N 
qubits in a translationally invariant state and determined 
the maximum possible nearest neigbour entanglement to 
be at least C = 0.434467 or t = 0.18876157 per qubit. 

The distributions for T3 and tab^ Fig- || (inset), were 
calculated numerically by samping a million pure three 
qubit states drawn randomly over the Haar measure, the 
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FIG. 2; Distribution of the 2-tangle for two (*), three (filled 
A) and four (■) qubit pure states; inset: distribution of the 
3-tangle (A) and pairwise tangle (filled A) for three qubit 
pure states, all sampled uniformly over the Haar measure. 



generalization of Eq. (||). The 2-tangle is now much 
more concentrated near zero than for two qubits, while 
the 3-tangle is broadly peaked around the average value. 

These values for pairwise and three-way tangles would 
be more useful if it was possible to continue in the same 
manner for four and more qubits. Currently, there is no 
known analytical expression for the 3-tangle of a mixed 
state of three qubits, nor is it known whether expres- 
sions equivalent to Eq. (^) can be found beyond three 
qubits. However, we can at least continue to evaluate 
pairwise entanglement for larger pure states. For four 
qubits the average pairwise entanglement evaluated nu- 
merically from a million randomly sampled four qubit 
pure states is 0.03138 ± 0.00006 tangle units per pair. 
Despite the low value for pairwise entanglement, a typ- 
ical four qubit pure state is still highly entangled over- 
all |23), as can be shown by considering the entropy of 
its subsystems. For a pure state of N qubits, the en- 
tropy Sn = 0. If the system is split into two pieces, of 
k qubits and {N — k) qubits, then each subsystem has 
the same entropy, Sk — 'S'(jv-fe) = X^i log2 A^, where 
the \i are the eigenvalues of the reduced density matrix 
Pk obtained by tracing out the other {N — k) qubits (or 
vice versa) . The entropy Sk measures how entangled the 
two subsystems are with each other. For four qubit pure 
states, {Si) = {S3) = 0.8661 ... and (^2) = 0.6653 . . . 
normalised per qubit such that < 5^ < 1. 

The four qubit pairwise distribution, shown in Fig. ||, 
is now strongly peaked around zero entanglement. A 
closer look at the numerical data reveals a new feature 
not present in the distributions for two or three qubit 
pure states: there is now a significant fraction of the pairs 
(24%) with zero entanglement, causing the outlying point 
at 67 for the bin that includes zero. This is not just zero 
to numerical accuracy. The formula for the concurrence 
of a pair of qubits in a mixed state p is given by Eq. 
(|l|). Clearly, if the values of the Xi are such that the 
maximum is zero by a significant margin, then the result 
is effectively exact. This trend continues for random pure 
states of five and six qubits, with (evaluated numerically 
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FIG. 3: Average negativity for subsets of m qubits selected 
from random pure states sampled uniformly over the Haar 
measure. The symbol gives m. For subsets of two qubits, the 
tangle (*) and concurrence (+) are also shown. 
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FIG. 4: Percentage of subsets of m qubits with no (free) en- 
tanglement in random pure states sampled uniformly over the 
Haar measure, where m is the symbol. 



over samples of 100,000) 80% and 99% respectively of 
pairs having zero entanglement in such states selected 
randomly over the Haar measure. This is the probably 
P2 of finding a chosen pair has zero entanglement. The 
probability Ps that the state as a whole has zero pairwise 
entanglement in all possible pairs is smaller, but also fast 
approaching 100%; by A'^ = 7 numerical sampling gives 
Ps > 99.9%. So we have a picture of typical pure states 
having less and less pairwise entanglement as the number 
of qubits increases, with the measure of states with some 
pairwise entanglement becoming essentially zero beyond 
7 qubits. In other words, the fall off is not gradual, but 
occurs sharply between 3 < iV < 7. 

This immediately begs the question, what happens to 
the entanglement in subsets of three or more qubits as 
N increases? While it isn't known how to calculate the 
3-tangle for mixed states of three qubits, there is another 
measure of entanglement that can be used to answer this 
question. The partial transpose of a density matrix ex- 
pressed in the standard basis is obtained by interchang- 
ing terms with a selected qubit in the opposite state Q . 
Such a matrix, denoted , can have negative eigenvalues 
if the original density matrix p is entangled. If has 
only positive eigenvalues for any combination of trans- 
posed qubits (up to half the total in p), then p is said to 
be PPT (positive partial transpose), and in most cases, 
p is separable. The exceptions have bound entanglement 
p5| , p6[ , which only occurs in systems with Hilbert space 
larger than 2 2. Bound entangled states are known to 
be relatively rare , of finite measure but smaller than 
the measure of separable states, which decreases expo- 
nentially in the number of qubits. 

A PPT test is sufficient to show that a three qubit sub- 
set has no free entanglement, but there is also a quanti- 
tative measure of entanglement that can be derived from 
p^ , the negativity. Following Zyczkowski ||l2| we define 
the negativity as twice the absolute value of the sum of 
the negative eigenvalues of p^; Vidal and Werner |2^ 
recently proved it is an entanglement monotone. Figure 
§ shows how the average values compare with the tan- 
gle and concurrence, two curves are shown for subsets 



of four and five qubits because there are different val- 
ues for transposing one qubit, or two qubits at once. A 
log-linear plot has been used to show that all the entan- 
glement measures fall faster than exponentially as the 
number of qubits in the pure state increases. 

Note, however, that the negativity does not satisfy any 
convenient quantitative relation like Eq. (^) for the tan- 
gle. It is thus not possible to identify finite amounts of 
three-way entanglement unless you also know that there 
is zero pairwise entanglement in the pure state. Con- 
versely, if the negativity of all subsets of three qubits 
is zero, there is still the possibility of bound entangle- 
ment. Nonetheless, since bound entanglement is a rela- 
tively rare phenomenon, if we are willing to neglect it in 
order to obtain the shape of the larger picture over all 
pure states, we can use the negativity to identify states 
that have approximately zero entanglement in subsets of 
three, four and five qubits. Fig. ^. In each case, it drops 
sharply to zero between N ~ 2m and N = 2m + 3, based 
on numerical analysis of random pure states sampled uni- 
formly over the Haar measure for TV < 13 and to < 6. 

This result can be explained at least qualitatively as 
follows. Page conjectured (later proved ^) that 
the average entropy of such a subsystem is (in the nota- 
tion of this paper) 

^kH2)-^J^ (5) 

where the approximate version is valid for 1 ^ fc < A^/2. 
These average entropies can also be calculated numeri- 
cally in the same manner as we calculated the tangle and 
concurrence, by sampling random pure states. Compari- 
son of these values with the theoretical formula provides 
a useful test of the accuracy of the numerical method, 
and agreement was generally better than 3 significant 
digits, even for relatively small sample sizes. A useful 
way of presenting {Sk) is shown in Fig. |^, where it is 
plotted per qubit for both the larger and smaller of the 
two subsystems. Equation (||) implies that the entropy of 
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FIG. 5: Average entropy of subsystems plotted per qubit for 
random pure states sampled uniformly over the Haar measure. 
Symbols give subsystem size, k and (A'^ — k). For N = 2, the 
average entropy corresponds to {Ef) for two qubits, as shown. 



the smaller subsystem is nearly maximal, which in turn 
implies that the subsystem is nearly maximally mixed, 
i.e. not entangled. In Fig. |5| this is shown by the lines 
for the smaller subsystem asymptoting to 1 as the dif- 
ference between k and TV becomes larger. There is a set 
of unentangled states of finite measure p2[ surrounding 
the maximally mixed state (p = 1). As N increases and 
the sirraller subsystem becomes more mixed, it enters this 
unentangled region for finite N, thus explaining the rapid 
transition to zero entanglement shown in Fig. |^. 

In conclusion, we have shown that in randomly cho- 
sen pure states of N qubits, there is a cut-off, m, be- 
low which subsets of qubits taken from the pure state by 
partial tracing over the remainder can be expected to be 
PPT (i.e. no free entanglement). For N < 13, the limited 
range accessible to our numerical studies, m > (iV — 3)/2. 
There are already experimental situations (such as ion 
traps) with more than six distinguishable quantum par- 
ticles in which such entangled states might be generated. 
Only looking for pairwise entanglement in such systems 
may give misleading results. 
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